Microwave response of thin niobium films under perpendicular static magnetic fields 



O 
O 

(N 



(N 



S3 
O 
O 

di- 
ss 

o3 



D. Janjusevic, M. S. Grbic, M. PozekQ A. Dulcic, and D. Paar 
Department of Physics, Faculty of Science, University of Zagreb, P. 0. Box 331, HR-10002 Zagreb, Croatia 

B. NebendahQ 

2. Physikalisches Institut, Universitat Stuttgart, D-70550 Stuttgart, Germany 

T. Wagner 

Max-Planck-Institute for Metals Research, Heisenbergstrasse 3, D-70156 Stuttgart, Germany 

The microwave response of high quality niobium films in a perpendicular static magnetic field 
has been investigated. The complex frequency shift was measured up to the upper critical fields. 
The data have been analyzed by the effective conductivity model for the type-fl superconductors 
in the mixed state. This model is found to yield consistent results for the coherence lengths in 
high-K superconducting samples, and can be used with HTSC even at temperatures much below 
T c . It is shown that for samples with high values of depinning frequency, one should measure both 
components of the complex frequency shift in order to determine the flow resistivity. 

The thick Nb film (160 nm) has low resistivity at 10 K, comparable to the best single crystals, 
and low k value. In contrast, the thinnest (10 nm) film has k « 9.5 and exhibits a high depinning 
frequency (~ 20GHz). The upper critical field determined from microwave measurements is related 
to the radius of nonoverlaping vortices, and appears to be larger than the one determined by the 
transition to the normal state. 

PACS numbers: 74.78.Db 74.25. Nf 74.25.Op 74.25.Qt 
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I. INTRODUCTION 

Microwave investigations of classical superconductors 
half a century ago have proven to be a very useful tool in 
the determination of intrinsic parameters in the Meiss- 
ner stat o 1 ^ 3 as well as the flow resistivity in the mixed 
statei 4 i 5 i 6 These investigations were also of large technical 
interest for the construction of microwave transmission 
lines and resonant structures. In the last twenty years 
high temperature superconductors (HTSC) have been ex- 
tensively investigated by microwave techniques. Among 
the greatest successes, one should mention the determi- 
nation of the temperature dependence of the penetration 
depth which evidenced the existence of nodes in the su- 
perconducting gap, leading to the d-wave explanation of 
HTSCii 

Flux flow resistivity can be determined from dc mea- 
surements when the current density J exceeds the crit- 
ical current density J c for vortex depinning. 8 However, 
in many cases of low and high temperature superconduc- 
tors, the dynamics of vortex motion is more complicated 
and cannot be interpreted just in terms of pinned and dc- 
pinned regimes. Thermally activated flux hopping may 
yield ohmic behavior even for current densities below J c . 9 
However, the resistivity observed in this regime should 
not be identified with that of flux flow. By increasing the 
current density there occurs a nonlinear transition to the 
flux flow regime. At much higher current densities one 
may encounter nonlinear effects in the flux flow leading 
ultimately to the point of instability where a dramatic in- 
crease of the measured voltage occursiiSiiiiiSii 3 ^^^ The 
nonlinear behavior is temperature and magnetic field de- 



pendent. Hence, the linear flux flow regime may occur 
inbetween the two nonlinear regimes, and may not always 
be well resolved. 

In contrast to dc resistivity, microwave measurements 
can be carried out with current densities much smaller 
than J c , and still yield flux flow resistivity. Namely, the 
vortices are not driven over the pinning potential barrier, 
but just oscillate within the potential well. 

Depinning frequency u>o separates two regimes of vor- 
tex oscillation. If the driving frequency lj is much larger 
than wo, the viscous drag force dominates over the restor- 
ing pinning force in the response of vortices to the mi- 
crowave Lorentz force. This is often the case with classi- 
cal superconductors which have ujq in the MHz range. In 
that case, the flux flow resistivity can be extracted from 
the microwave absorption curves only. When loq w lj, 
or higher, the restoring pinning force is comparable to 
the viscous drag force, and the response of the vortices 
is more complex. This is typical of HTSC where the 
depinning frequency is found to be of the order of 10 
GHz, 1 ?" but may be found also in low temperature super- 
conductors such as in very thin Nb film analyzed in this 
paper below. In such cases, one needs both, microwave 
absorption and dispersion to determine the depinning fre- 
quency and flux flow parameters. The analysis usually 
employs the models of effective conductivity in the mixed 
stat o 17 i 18 i 19 based on the Bardeen-Stephen models 

The real values of B C 2 in HTSC remain experimentally 
unreachable except for the narrow temperature range be- 
low T c , and the B c i values extracted from the effective 
conductivity data could not be experimentally verified. 
In order to probe the effective conductivity model, we 
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have performed a series of measurements on thin films 
and a single crystal of niobium, a classical type-II su- 
perconductor. The comparison of thick niobium samples 
with HTSC samples is not ideal since the upper criti- 
cal fields and depinning frequencies in niobium are much 
lower and the vortex distance doesn't always allow the 
use of effective conductivity models. However, as the 
film thickness is reduced, one expects considerable en- 
hancement of B C 2, due to the reduced coherence length, 
and much stronger pinning, due to the surface effects. 

We have measured high quality niobium thin films us- 
ing a cavity perturbation method. The temperature and 
field dependence of the complex frequency shift have been 
measured up to the upper critical fields in order to test 
the validity of the effective conductivity models. The re- 
sults of the present analysis should be taken into account 
when the mixed state of HTSC is investigated by the 
microwave methods. 



II. SAMPLES 
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FIG. 1: (Color online) Thickness dependence of the transi- 
tion temperature of the measured niobium films (full squares). 
The data are compared with clean samples of Wolf et al»— 
and Kodama et ali-^ (dashed line), and with dirty samples of 
Minhaj et all— (circles). 



The Nb films were deposited via molecular beam epi- 
taxy (MBE) in a commercial system from DCA instru- 
ments (Finland). The base pressure of the system is 10 -9 
Pa. The (0001) surfaces of the sapphire substrates (a- 
AI2O3) were prepared by sputter cleaning with Ar ions (1 
keV) and subsequent annealing at 1000°C in ultra-high 
vacuum (UHV). This annealing procedure is necessary to 
remove the embedded Ar gas atoms and to recover the 
sapphire surface crystallography and morphology. De- 
tails of the surface preparation procedure can be found 
elsewhere i^iiSSiS The substrate treatment resulted in un- 
reconstructed surfaces, as revealed by reflection high en- 
ergy electron diffraction (RHEED). Niobium (4 N pu- 
rity) was evaporated from an electron beam evaporator 
(substrate temperature 900°C; film thicknesses 10 nm, 
40 nm, 160 nm) and the growth rate was monitored using 
a quartz oscillator. Typical growth rates were between 
0.01 and 0.05 nm/s. As revealed by in-situ RHEED inves- 
tigations, the Nb films grew epitaxially on the (0001) sap- 
phire substrate^ A 2 nm thick protective layer of SiC>2 
has been deposited onto each film in order to prevent 
niobium oxidation (physical vapour deposition at room 
temperature). Transmission electron microscopy (TEM) 
investigations have shown that the film of nominal thick- 
ness 40 nm is actually 36 nm thick, indicating that the 
other films are also slightly thinner than their nominal 
thickness. The samples were cut by a diamond saw into 
pieces 3x0.5 mm 2 (40 nm and 160 nm films), and 2.3x1 
mm 2 (10 nm film), suitable for cavity perturbation mea- 
surements. High purity single crystal of niobium was pur- 
chased from Metal Crystals & Oxides Ltd., Cambridge. 
Its dimensions were 3x2x0.5 mm 3 . 

Superconducting properties of thin niobium films 
strongly depend on their thickness, purity and prepa- 
ration conditionsiSiiSSSiSL 2 ^ 9 . In order to compare our 
samples with those measured earlier ^2^^ W e have de- 



termined T c of each film by the so called 90% criterion, 
i.e. the transition temperature was defined as the tem- 
perature where microwave absorption A(1/2Q) reaches 
90% of its normal state value. For the samples of the 
same composition and purity one expects that T c is re- 
duced aproximately with the inverse thickness d _1 due to 
the proximity effects^ weak localization and increased 
residual resistivity* 2 ^ We show in Fig^ the comparison 
of transition temperatures of our thin films with results 
of some other authors. The T c values of present films are 
similar to the films measured by Gubin et al.^ and, if 
extrapolated to ultrathin 2 nm, with the film measured 
by Hsu and Kapitulnik. 29 . Obviously, the films measured 
in this paper are of good quality. 

To extract various superconducting parameters one 
needs three quantities: T c , p n and S = —dB C 2/dT at 
T c . The B C 2 values determined also by the 90% crite- 
rion from our measurements gave the S'-values 0.47 T/K, 
0.11 T/K and 0.08 T/K for 10 nm, 40 nm and 160 nm 
samples, respectively. 

The resistivities of the films have been determined by 
standard four contact method. Their residual values at 
T = 10 K were p„(10 nm) = 15.2 ^Slcm, p„(40 nm) = 
1.4 /iOcm and p n (l60 nm) = 0.33 /if2cm, respectively. 
The residual resistivity of 160 nm film was found to be 
an order of magnitude lower than in the film studied by 
Gubin et al.j2£ and comparable to the best bulk single 
crystals showing that the films grew with perfect epi- 
taxy with virtually no defects. The residual resistivities 
of thinner samples are, therefore, increased only due to 
the surface scattering, but the values are still below the 
values reported by Gubin et ali2& For the 10 nm film we 
estimate k 9.5, for the 40 nm film k « 1.5, and for 
the 160 nm film we can take the single crystal value of 
K« 0.9. 
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III. EXPERIMENTAL DETAILS 

Microwave measurements were carried out in a high- 
Q elliptical cavity made of copper resonating in e TE m 
mode at 9.3 GHz or in e TE 113 mode at 17.5 GHz. For 
both modes used, microwave magnetic field has node in 
the center of the cavity while microwave electric field is at 
maximum. The sample was mounted on a sapphire sam- 
ple holder and placed in the center of the cavity where 
the microwave electric field E mw has its maximum in 
both modes. The sample was oriented with the longest 
side parallel to E rnw . Experimental setup included Ox- 
ford Systems superconducting magnet with B^c -L E rnw . 
While E mw is always in the film plane, B c i c can be in 
the film plane or perpendicular to it. Directly mea- 
sured quantities are the Q-factor and the resonant fre- 
quency / of the cavity loaded with the sample. The Q- 
factor was measured by a modulation technique described 
elsewhere^ The empty cavity absorption (1/2Q) was 
substracted from the measured data and the presented 
experimental curves are due to the samples themselves. 
An automatic frequency control (AFC) system was used 
to set the source frequency always in resonance with the 
cavity. Thus, the frequency shift can be measured as the 
temperature of the sample or static magnetic field is var- 
ied. The two measured quantities represent the complex 
frequency shift Au/uj = Af/f + zA(l/2Q). 



conductivityiii£ii£ which is a combination of normal 
conductivity in the vortex cores and the conductivity of 
the condensed electrons outside the cores. The effective 
conductivity in an oscillating electric field is given by: 



<7 e ff (1 - &)(oi - ia 2 ) + ba n u n 1 - i(u) /u) 

The first term is due to the microwave current outside the 
vortex cores, and the second is due to the normal current 
in the cores of the oscillating vortices. The meaning of 
the parameter b in Eq. is the volume fraction of the 
sample taken by the normal vortex cores. This parameter 
determines the resistivity in the flux flow regime p/ / p n ~ 
The depinning frequency u>q may vary, depending on 
sample, field and temperature from strongly pinned case 
(luq ^> v) to the flux flow limit (u>q <C oS). In Eq. (0) 
the zero field conductivity is a\ — ia 2 , and a n is the nor- 
mal state conductivity. The model is limitted to high 
K-values and to magnetic fields much lower than upper 
critical fields, where vortices don't overlap significantly. 

From the experimentally obtained field dependent 
complex conductivity extracted using Eqs. (TJ) and © 
one can, by numerical inversion of Eq. , determine the 
values of b and luq/lo. 



IV. EFFECTIVE CONDUCTIVITY 

For the extraction of complex conductivity data from 
the measured complex frequency shift of a thin super- 
conducting sample in the microwave electric field, one 
can utilize the general solution for the complex frequency 
shift by Peligrad et The shift from a perfect con- 

ductor state is given by 
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tanh(ifcd/2) 
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where d is the thickness of the slab, TV is the depolariza- 
tion factor and T is the dimensionless filling factor of the 
sample in the cavity. The complex wave vector k is given 
by 



k = k i jur 



-'-) 



(2) 



where ko = u^/Ji^e® is the vacuum wave vector. It de- 
scribes generally any set of material parameters. For a 
nonmagnetic metal one can take JT r — e r = 1 and the 
main contribution to k comes from conductivity. Using 
these equations one can by numerical inversion of experi- 
mentally obtained complex frequency shift determine the 
complex conductivity a of the sample. 

For a superconducting sample in the mixed state, in 
intermediate fields one can define the effective complex 



V. RESULTS AND DISCUSSION 

The measured complex frequency shift of the three 
films and the single crystal is presented in Fig. [21 for 
various values of DC magnetic field and for a driving 
frequency 9.3 GHz. One may notice a huge difference 
in signal intensities between thin and thick samples due 
to very different filling factors T. Namely, a conducting 
sample profoundly changes the electric field at the centre 
of the cavity with respect to the empty cavity field. The 
sample acts as a partial short for the electric field lines. 35 
The field outside of a conductor is stronger at rear and 
front sides, and this enhancement is larger for thinner 
films. As a consequence, the microwave measurements 
in electric field are very sensitive for thin films, while for 
thicker samples the signal/noise ratio is reduced. 

Looking at the zero field measurements in Fig. [21 
one observes a smooth transition to the superconduct- 
ing state for the thinnest film, and sharp transitions for 
thicker samples. In the 10 nm sample the fluctuation ef- 
fects become considerable since the thickness becomes 
lower than the coherence length. 

From the measured frequency shift and known nor- 
mal state resistivity p n one can determine the complex 
conductivity in the whole temperature range and con- 
sequently the zero temperature penetration depth. We 
have obtained A(T = 0) = 285 nm for the 10 nm film, 
comparable to the result of Gubin et ali2& for the film 
of the same thickness. The filling factor T was deter- 
mined at T — 10 K from the measured frequency shift 
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FIG. 3: (Color online) Field dependence of the complex fre- 
quency shift for the 10 nm film at various temperatures. 
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FIG. 4: (Color online) Parameters b and ujo/lu of the effective 
conductivity model extracted from the data of Fig. |^1 The 
dashed line shows the extrapolation of the linear section at 
T = 2.9K. 



10 nm 
T=2.9K 




0.0 0.1 



t — 1 — i — ' — r 

0.2 0.3 0.4 0.5 0.6 



S(T) 



FIG. 5: (Color online) Comparison of b obtained by full in- 
version of complex frequency shift (full black circles) and b 
obtained from absorption measurement alone under assump- 
tion that depinning frequency is much lower than driving fre- 
quency (red open circles): a) 10 nm film at 2.9 K; b) 40 nm 
film at 5 K. 



this linear section, shown by the dashed line in Fig. 0Ji, 
would give B* 2 (2.9 K) = 1.9 T, and £(2.9 K) = 41 nm. 
At higher fields b changes its slope and approaches the 
normal state value at field B C 2 lower than B* 2 . It is the 
region where vortices start to overlap and the supercon- 
ducting order parameter is reduced throughout the sam- 
ple. The deviation from linearity starts roughly at the 
field B — 0.63 B c2 , where the distance between vortex 
centers is 3.4 times the coherence length. 

One should mention that numerically determined b val- 
ues don't always reach unity at B C 2. Close to B C 2, b looses 
its meaning as a parameter of the effective conductivity 
model. There are no more well defined vortices and the 
conductivity is dominated by fluctuations. One can still 
try to numerically obtain b and u>q, but this leads to neg- 
ative values of depinning frequencies, clearly showing the 
failure of the effective conductivity model in that region. 

The determination of flow resistivity would not be pos- 
sible if only the microwave absorption had been mea- 
sured. The effective conductivity depends on two param- 
eters (b and wo), and one has to measure two quantities, 
absorption and real frequency shift, in order to deter- 
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FIG. 6: (Color online) Field dependence of the depinning fre- 
quency vq — u>o/2n obtained by driving frequency 9.3 GHz 
for three niobium films at T/T c ~ 0.5: 10 nm film at 2.9 K 
(black full circles) ; 40 nm film at 5 K (red opened triangles) 
and 160 nm film at 5 K (green full diamonds). The depin- 
ning frequency obtained by driving frequency 17.5 GHz for 
the 10 nm film at 2.9 K is shown by open black circles. 



mine both parameters. If one supposed that the pinning 
at microwave frequencies were negligible, one could have 
determined the flow resistivity from absorption measure- 
ments alone. This would lead to wrong conclusions. We 
show in Fig. |U< the comparison of the two ways of rea- 
soning for the 10 nm film at 2.9 K. The full circles show 
the parameter b determined from the complex frequency 
shift with the depinning frequency as the second result of 
the numerical inversion. The empty circles are the result 
of the numerical inversion of the microwave absorption 
only, under the assumption that loq = 0. One can see 
that the results differ strongly in the low field region, i.e. 
in the region where the pinning is strong. When the de- 
pinning frequency falls well below the driving frequency, 
the two curves overlap. A similar analysis for the 40 nm 
sample is shown in Fig. [SJ}. The two ways of reasoning 
give closer results since the depinning frequency for this 
sample is lower than the driving frequency (u>o < 0.7u>, 
see Fig. |SJ . 

If one is dealing with the samples whose upper critical 
field is experimentally unreachable (as in HTSC), one is 
lead to infer the upper critical field from the slope of 
b in low fields. Taking the flow resistivity values from 
microwave absorption, without taking the real frequency 
shift into account, one could reach erroneous values of 
the upper critical field and GL coherence length. This is 
especially true for samples with high values of depinning 
frequency, which is usually the case in HTSC. 

Bulk metallic superconductors usually don't show pin- 
ning at microwave frequencies. 6 In Fig. we show the 
depinning frequency of our three niobium films at T « 
0.5 T c . The depinning frequency for the 10 nm film at 
low fields is slightly above 20 GHz, and for the 40 nm 
film it is approximately 7 GHz. For the 160 nm film it 
is below 1 GHz, i.e. in the frequency range that would 
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FIG. 7: (Color online) Upper critical field B C 2 of the 10 nm 
film obtained by the 90 % criterion from the microwave ab- 
sorption measurements (black full circles) . Upper critical field 
B* 2 of the same sample obtained from the linear section of b 
(black opened circles) . The same is shown for the 40 nm film 
by red triangles. 



be better examined by lower driving frequency. Obvi- 
ously, in thin films placed in a perpendicular magnetic 
field, vortex pinning is dominated by surface pinning cen- 
ters. Thicker films approach bulk values of depinning 
frequencies, showing that the presently grown niobium 
films don't have bulk pinning centers, which are absent 
also in the best single crystals. 

To check the consistency of depinning frequency values 
obtained by different driving frequencies, we obtained the 
depinning frequency of 10 nm sample from measurements 
in the e TEn3 mode with driving frequency of 17.5 GHz. 
The result is shown by the open circles in Fig. [fjl with 
satisfactory agreement between the two measurement se- 
ries. 

The above analysis for the 10 nm film clearly shows 
linear section of b in moderate fields. Already for the 
40 nm film, the linearity is not perfect. For the 160 nm 
film and for single crystal there were no linear sections. It 
is not surprising since the simple Bardeen-Stephen model 
predicts a linear dependence of Pf/p n on B only in the 
high-re limit. Obviously, only the thinnest film considered 
here fulfills the necessary conditions for the analysis in 
terms of the effective conductivity model. This analysis 
is a good model for the analysis of high-re superconductor 
with high values of upper critical fields, typically HTSC. 

One can compare the values of B c2 , obtained by the 
90 % criterion, with the B* 2 values obtained from the 
linear section of field dependence of b by taking into ac- 
count the Larkin-Ovchinnikov correction b = 0.9 B/B* 2 . 
The comparison is shown in Fig.0 Full circles show B c2 
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obtained by the 90 % criterion and opened circles show 
B* 2 values obtained from linear section of b. One can see 
that B* 2 is always slightly higher than B C 2 which would 
be measured by direct transport or magnetization mea- 
surement of the crossover to the normal state, if exper- 
imentally reachable. The reduction of order parameter 
due to the overlapping of vortices in higher fields is one 
of the reasons. The other reason for further suppression 
of the actual B C 2 value can be found in spin paramagnetic 
effects and spin-orbit scattering^ This effect is not sig- 
nificant in our niobium samples since the upper critical 
fields are low enough, but in HTSC, one can expect sig- 
nificant suppression of B C 2 in the fields of the order of 
100 T. As a consequence, the transition to the normal 
state is due to a combination of effects, and could not be 
related uniquely to the coherence length. However, the 
microwave determination of flow resistance in low field 
region gives the best estimate of the coherence length. 

The comparison of B C 2 and B* 2 for 40 nm film is also 
shown in Fig. [7| The B* 2 values are significantly higher 
than B C 2, but one should take this result with caution 
since there is no strictly linear part of b field dependence. 
This is not surprising since k ~ 1.5 doesn't meet a high-K 
condition needed for B* 2 extraction. 

VI. CONCLUSIONS 

We have measured temperature and magnetic field 
dependence of the microwave response of high quality 
niobium films in a perpendicular static magnetic field. 
The data have been analyzed by the effective conductiv- 
ity model for the high-K type-II superconductors in the 
mixed state. Our analysis shows that this model yields 
consistent results when applied in fields not too close to 
the transition to the normal state. Thin niobium films 
show strong pinning probably due to the surface pinning 
centers. The depinning frequency for the 10 nm film is 
ujq ss 20 GHz, and for the 40 nm film it is loq « 7 GHz. 
In thicker samples the depinning frequency is lower than 
1 GHz. One should measure both, the microwave ab- 
sorption and the fequency shift in order to determine 
flow resistivity for samples with high values of depinning 
frequency, as usually found in HTSC. 

Although the bulk niobium samples have GL parame- 



ter k of the order of unity, very thin niobium films have 
higher K-values. They can serve as model system for the 
upper critical field determination of high-K superconduct- 
ing samples. In HTSC generally, upper critical fields are 
of the order of hundreds tesla, and they are unreachable 
by laboratory magnets. Conventional methods for the 
determination of B C 2 (transport measurements, suscepti- 
bility, torque,...) are based on the crossover to the normal 
state when the applied field equals to B c2 . The determi- 
nation of upper critical field is, therefore, reduced to a 
narrow temperature region close to T c . The values for 
lower temperatures are usually determined by an extrap- 
olation. The correctness of such extrapolations depends 
on the model used and there is no assurance that it holds 
also in HTSC. On the other hand, pulsed methods used 
to reach high applied fields (up to 400 T— ) are nonequi- 
librium methods and very sensitive to induction, criti- 
cal currents etc. In this paper microwave response was 
measured up to the upper critical field which is for nio- 
bium samples reachable by conventional magnets. The 
upper critical fields are determined by the flow resistiv- 
ity in low fields {B* 2 ) and by the crossover to the normal 
state (B C 2). The two methods give slightly different re- 
sults. The B C 2 defined as the crossover field to the nor- 
mal state has certainly great technological importance. 
However, at high fields, typical for HTSC, paramagnetic 
effects induce depairing and reduce order parameter. If 
one calculates the coherence lengths from B C 2 values ob- 
tained in this way, the values are overestimated. For the 
considerations of the physical parameters of the super- 
conducting state, it is better to determine the coherence 
length in rather low magnetic fields, deeply in the su- 
perconducting state. The microwave method tested here 
is certainly suitable for the determination of B* 2 which 
gives correct values of effective vortex radii at a given 
temperature in the fields low enough that vortices don't 
overlap. Therefore, the method is relevant for the deter- 
mination of coherence lengths, especially in HTSC. 
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